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SUMMARY 

This paper presents a comprehensive approach in analyzing the dynamic 
behavior of multistage gear transmission systems with the effects of gearbox 
induced vibrations and mass imbalances of the rotor. The modal method, with 
undamped frequencies and planar mode shapes, is used to reduce the degrees of 
freedom of the gear system for time-transient dynamic analysis. Both the 
lateral and torsional vibration modes of each rotor-bearing-gear stage as well 
as the interstage vibrational characteristics are coupled together through 
localized gear mesh tooth interactions. In addition, gearbox vibrations are 
also coupled to the rotor-bearing-gear system dynamics through bearing support 
forces between the rotor and the gearbox. Transient and steady state dynamics 
of lateral and torsional vibrations of the geared system are examined in both 
time and frequency domains to develop interpretations of the overall modal 
dynamic characteristics under various operating conditions. A typical three- 
stage geared system is used as an example. Effects of mass imbalance and 
gearbox vibrations on the system dynamic behavior are presented in terms of 
modal excitation functions for both lateral and torsional vibrations. Opera- 
tional characteristics and conclusions are drawn from the results presented. 


INTRODUCTION 

In the continuous search for the improvements in operational life, effi- 
ciency, maintainability, and reliability in gear transmission systems, one of 
the major objectives in design improvement is the reduction of noise and vibra- 
tion in the transmission system. Two main streams of work have been carried 
out in the areas of noise and vibration reduction; namely (1) the localized 
gear-tooth stress/thermal effects during gear interactions [Cornell 1981, Lin 
1988, Boyd 1987, Savage 1986], and (2) the overall global dynamic behavior 
[August 1986, Choy 1989, David 1987, 1988, Mitchell 1985] of the transmission 
systems . 

The work presented in this paper is the development and application of a 
comprehensive approach in simulating the overall dynamics of a multistage rotor- 
bearing-gear system. The modal method is applied to reduce the degrees of 
freedom of the global system of equations. Various operational characteristics 
are simulated in the analysis, such as 



(1) Multistage rotor-gear configuration 

(2) The coupling of gearbox vibrations 

(3) The effects of rotor mass- imbalance and shaft bow 

(4) The effects of inertia-gyroscopic moments 

(5) The effects of nonlinear gear mesh stiffnesses 

(6) Variable shaft geometry and bearing support characteristics 

The modal equations of motion are solved to evaluate system acceleration at 
each time step. A self-adaptive variable time-stepping integration technique 
[Choy 1987, 1988, 1989] is used to calculate the transient dynamics of the 
system. A typical three-stage rotor-bearing-gear system is used as an example 
in this analysis. Results are presented in both time and frequency domains to 
develop a vibration signature analysis of the system. 


NOMENCLATURE 


Ai(t) 

modal function of the i*“ mode in x-direction 

Ati(t) 

modal function of the i*h mode in 0-direction 

B i ( t ) 

modal function of the i*h mode in y-direction 

[C T ] 

torsional damping matrix 

» I^xy] 

bearing direct and cross-coupling damping matrices 

[CyxJ f [Cyy] 


FCttO 

gear mesh torque 

F Gx( t ) > F Gy( 4 ) 

gear mesh force in x- and y-directions 

F t (t) 

external excitation moment 

Fx( 1 ) >Fy( t ) 

external excitation forces 

[ 3 ’ f Gy ] 

gyroscopic-inertia matrix 

[I] 

identi ty matrix 

[J] 

rotational mass moment of inertia matrix 

K s 

shaft stiffness matrix 

[K T ] 

torsional stiffness matrix 

Ktki 

gear mesh stiffness between i*h and Ic** 1 rotor 

[ ^xx 3 > t ^xy 3 

bearing direct and cross-coupling damping matrix 

[ Kyx 1 > 1 Kyy 1 


[M] 

mass-inertia matrix 


2 


^ci 

tf 

Xci.Yd 

X F ,Y F 

a ki 

[A 2 ], [A 2 ] 
[*lk*[*tlk 


radius of gear in the i** 1 rotor 
gear generated torque 

gear displacements in x- and y-directions of the i** 1 rotor 

gear forces in x- and y-directions 

angle of tooth mesh between k** 1 and i^ rotor 

lateral and torsional eigenvalue diagonal matrices 

lateral and torsional orthonormal eigenvalue diagonal matrices 
of the k*h rotor 


DEVELOPMENT OF EQUATIONS OF MOTION 

The equations of motion for a single stage multimass rotor-bearing-gear 
system with the effects of gearbox motion induced vibrations at the bearing sup- 
ports, rotor inertia and gyroscopic effects, and excitations from rotor mass 
imbalance and residual runouts can be written in matrix form [Choy 1987, 1989] 
for the X-Y plane as 

[M]{X> * [G x ](X| ♦ [C^HX - X b } + (C xy ]{Y - Y b } ♦ [G y ]{Y) ♦ [K^ ♦ KJfX} 

- - [X^HY - Y b } - {F x (t)J + fF Cx (t)) (1) 
and in the Y-Z plane as 

[M](Y} - [C X ]{Y} * [CyJlX - X b } * [C^HY - Y b ) - (G y ]{Y} ♦ [K yy * K s ]{Y} 

- [KyjfYh) ♦ [K^HY - Y b ) - (F y (t)J + (F 0y (t)| (2) 

Here F x and Fy are force excitations from the effects of mass imbalance, 
shaft residual bow and bearing support base motion in both X- and Y-directions. 
The X and Y gear mesh forces, Fq x and FQy, are induced from the gear teeth 
interaction with other coupled gear stages, Tne bearing forces are evaluated 
through the relative motion between the rotor {X}, {Y} and the gearbox {X^} , 

{Yt,} at the bearing locations [Choy 1987], The mass-inertia and gyroscopic 
effects are incorporated in the mass matrix [M] and the gyroscopic matrices 
[G x ] and [Gy], The coupled torsional equations of motion for the single rotor- 
bearing-gear system can be written as follows: 

[JH®} + [C T ]{0} + [1C,. ]{©> = {F T (t)> + (F Gt (t)} (3) 

In equation (3), (Fx(t)} represents the externally applied torque and {Fct(t)} 
represents the gear mesh induced moment. Note that equations (1) to (3) 
repeat for each single gear/rotor stage. The gear mesh forces couple the force 
equations of each stage to each other as well as the torsional equations to the 
lateral equations [Choy 1989, Cornell 1981, David 1987, 1988]. The coupling 
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relationships between the torsional and the lateral vibrations and the dynamics 
of each individual gear/rotor stage are derived in the next section. 


COUPLING IN GEAR MESHES 

The torsional and lateral vibrations of a single, individual rotor and the 
dynamic relationships between each gear stage are coupled through the nonlinear 
interactions in the gear mesh. Gear mesh forces and moments are evaluated as 
functions of relative motion and rotation between two meshing gears and the cor- 
responding gear mesh stiffnesses. These gear mesh stiffnesses vary in a repeat- 
ing nonlinear pattern [August 1986, Cornell 1981, Savage 1986] with each tooth 
pass engagement period and can be represented by a high order polynomial 
[Cornell 1981, Boyd 1987], The periodic nature of such nonlinear mesh stiff- 
nesses can also act as a source of steady state excitation to the gear system, 
With the coordinate system as shown in figure 1, the following gear mesh cou- 
pling equations can be established by equating force and moment [Choy 1989]. 

For the k stage gear of the system, summing force in the X-direction results 
in the following: 

n 

X F = X K tki [ ~ R ci 0 ci + R ck 0 ck + < X ci " Y ck )cosa ki 

i =1 , i*r 

+ (Y ci - Y ck* sina ki lcos<I ki (4 > 


Summing force in the Y-direction results in 
n 


y f 


E K tkiC- R ci 9 ck - R ck 9 ck + < X ci - X ck )cos “ki 
i=l , i*k 

* < Y ci - Y ck )si "“ki Isina ki 


(5) 


Summing moment in the Z-direction results in 

n 

T F = X R ck [K tki ( - R ci 0 ci ~ R ck 0 ck )] 
i=l , i*k 


where n is the number of stages in the system. 


MODAL ANALYSIS 

To reduce the computational effort, the number of degrees of freedom of 
the system of equations is reduced through modal transformation. Orthonormal 
modes for each individual rotor-bearing stage are obtained by solving the 
uncoupled system homogeneous characteristic equations. With the modal expansion 
approach [Choy 1987, 1988, 1989], the motion of the system can be expressed as 
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( 7 ) 


m 

W - A i{+i> 
i=l 

m 

{Y} = ^ Bi{<t>i> (8) 

i=l 

ro 

{0} = ]T At^^} (9) 

i = l 

where m is the number of modes used. The orthogonality conditions of the 
modes can be expressed as 

[4>] T [K][*] = [A 2 ] (10) 

[* t ] T [K t ][* t ] = t A tl (11) 

[<I>] T [M][^] - [* t ] T [J][<l> t ] = [I] (12) 

With the modal expansion and the orthogonality conditions, the modal equations 
of motion [Choy 1989] can be written as follows: 

X-Z equation 

{A} + [4>] T [G x ][<t>]{A} + [<j>] T [C xx ][4>]{A} + [<K] T [C xy ][<l>]{B} + [*] T [G y ] [4>] {B} 

+ (A 2 ]{A} + [<f] T [K xy ][4»]{B} = [<I>] T {F x (t) + F Gx (t)} (13) 

Y-Z equation 

{B} - [*] T [G x ][ 4>]{B} + [$] T [C yx ][<l»]{A} + [4>] T [C yy ][<l»]{B} - [<t>] T [G y ] [<J>] {A} 

+ [A 2 ]{B} + [<l>] T [K yx ][$]{A} = [4>] T {F y (t) + F Gy (t)} (14) 

0-equation 

{A t } + [4> t ] T [ C T ] [<J> t ] {A t } + [A 2 ]{A t } = [* t ] T {F t (t) + F ct (t)} (15) 

The gear mesh forces and moment can also be expressed in the modal form, for 
the k** 1 stage, as follows: 
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m 


- Z+kji 


j-i 


n 


X K tki [ - R ci e ci - R ck e ck + (x ci - X ck } 
i=l , Uk 


cosai . 
ki 


+ (y ci - y ck> sina ki ]c0Sa ki ( < 16 ) 


ID 


[ * ] k { V = S *kjl X K tki [ - R ci e ci - R ck e ck + (x ci - x ck )cosa ki 
j = l I i=l ,l-,k 

* (y ci - y ck* sin «ki I < 17 > 


(4 ' 1 k |T F> -Z* kiI E R ckf K tki(- R ci e ci - R ck 9 ck)l 
j=l i 1=1 , 1-k 


( 18 ) 


where K is the stage number, j is the mode number, and 1 is the station 
location of the gear mesh. 


SOLUTION PROCEDURE 

Rearrange the modal equations of motion developed in equations (13) 
to (15) into the following: 

X-equat ion 

{A} = -[*] T [G x ][<D]{A} - [$] T [C xx ][<t>]{A} - [<t>] T [C ] [$]{!)} - [<t>] T [C ][<t] {B} 

y y 

- [A 2 ] {A} - [<t>] T [K xy ][<t»]{B} + [4>] T {F x (t) + F Gx (t )} (19) 

Y-equat ion 

{B} = [<*] T [G x ][<t]{B} - [ < t ) ] T [C yx ] [4>] {A} - ['t>] T [C yy ][d>]{B} + [4*] T [G y ] [$] {B } 

- [A 2 ] {B} - [<t>] T [K yx ] [4.] {A} + [<f] T {F y (t) + F Gy ( t )} (20) 

©-equation 

{A t } = — [ < J» T ] T [ C T ] [<t> t ] { A t } - [A^]{A t }[4. t ] T + (F t (t) + F ct (t)} (21) 

A variable time stepping Newmark-Beta integration scheme evaluates the modal 
velocity and displacement at each time interval. In turn, equations (7) to (9) 
transform the modal displacements into absolute relative displacements in fixed 
coordinates. The gear mesh forces can be evaluated by the relative motion 
between the gear teeth using the nonlinear stiffnesses developed for the gear 
mesh interaction. The effects of gearbox vibration can be calculated as nonlin- 
ear bearing forces through the relative vibration between the gearbox and 
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the shaft at the bearing locations. Since the gearbox and the shaft both are 
vibrating independently, a separate transient integration scheme is required 
for each system at every time step before the coupling of bearing forces. 


DISCUSSION OF RESULTS 

To demonstrate the application of the discussed analytical method, a typi- 
cal three-stage gear transmission given in figure 2 is used as an example. All 
three gear stages have an identical 36-tooth gear and a mesh contact ratio of 
1.6. Stage I of the system is the driver with a rotational speed of 3000 rpm 
and an input torque of 2.25 kN-m. Both stages I and II are supported by three 
bearings while stage III only has two bearing supports. The first two bearings 
in stages I to III and the third bearings in stages I and II are identical. 

The rotor in stages I and II is 2.4 in. (6.1 cm) in diameter and has a length 
of 40 in. (101 cm) and 30 in. (76 cm), respectively. The rotor in stage III is 
smaller with a diameter of 1.4 in. (3.5 cm) and a length of 20 in. (51 cm). 

For demonstration purposes, the first three basic modes of the system are used 
in the modal analysis. Table I provides the information on the modal frequen- 
cies of the system in both lateral and torsional directions. Figures 3 and 4 
give the orthonormal mode shapes of the first two basic torsional and lateral 
modes for all three rotor stages. Note that the general form of mode shapes in 
stage I is more like a mirror image compacted to those of stages II and III 
which are generally identical. This is due to the fact that a nominal support 
stiffness is assumed at the gear locations, which produce a substantial differ- 
ence on the vibratory modes of each rotor stage. 

In order to investigate the effects of gearbox vibrations and rotor mass 
imbalance, results from four major cases of external excitations are examined 
in this study, namely 

Case 1 . - No gearbox vibration and small mass imbalance 
Case 2. - Gearbox vibration and small mass imbalance 
Case 3 . - No gearbox vibration and large mass imbalance 
Case 4. - Gearbox vibration and large mass imbalance 

The gearbox vibrational effects are assumed to be a steady state vibration 
motion of 1200 Hz at the bearing supports. Vibrational orbits of all three 
rotor stages at the gear location for all four cases, respectively, are given 
in figures 5 to 7. Note that the vibratory characteristics of stage I is more 
influenced by mass imbalance while those of stages II and III are affected by 
motion of the gearbox through the bearing supports. As we can see from fig- 
ure 5(b), the effects of bearing support motion to stage I gear vibrational 
orbits are quite small. This is due to the fact that the rotor of stage I is 
substantially more flexible than the other two rotors and less influence from 
the bearing can be effectively converted to the system. In addition, although 
similar mass imbalances are applied to each rotor stage, vibratory orbits in 
stages II and III are smaller in amplitude because of their larger rotor 
stiffnesses . 

Figures 8 to 10 depict the first modal component of rotor stages I to III, 
respectively. Figure 8 shows the effects of the four external excitations on 
the first modal vibratory component of stage I in both time and frequency 
domains. Note that figures 8(a) and (b) show that, for small imbalance with no 
gearbox motion, only very small vibratory motion will result except for the 
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dc component at zero frequency. Figures 8(c) and (d) represent the first mode 
rotor vibrations resulting from gearbox motion. Note that the major component 
occurs at 1200 Hz (input vibratory frequency at the bearing supports) with a 
minor component at 115 Hz (first natural frequency of stage I). Figures 8(e) 
and (f) show the effect of mass imbalance, which results in a large component 
at the rotor speed of 50 Hz and a smaller component at the first mode of 115 Hz, 

Figures 8(g) and (h) produce the combined effect of gearbox vibration and 

imbalance on rotor first mode vibrations. By examining the relative amplitudes 

of the major components, rotor speed, and gearbox vibration frequency, with 

their corresponding input excitation magnitude, the sensitivity of the rotor to 
imbalance and gearbox vibration can be determined. Similar conclusions can be 
arrived in figures 9 and 10 for stages II and III, except that the effects of 
mass imbalance on first mode vibration is significantly reduced. As we can see 
in figures 9(g) and (h) and 10(g) and (h), when both excitations are present, 
the gearbox motion has a much more dominating affect in rotor vibration for 
stages II and III as compared to stage I. Figures 11 to 13 show the vibration 
of the second mode component for the three rotor stages, respectively. Note in 
figure 11, for stage I the second mode vibrational characteristics are very sim- 
ilar to those of the first mode. However, magnitudes of the second mode are 
approximately 20 percent of those of the first mode. This is because the first 
mode in the first stage is more easily excited by the imbalances and gear forces 
acting at the gear location. Further examinations of figures 12 and 13 show 
that the second mode is more excited than the first mode in the second and 
third rotor stages. A more detailed study of figures 3 and 4 reveals that the 
second mode in stages II and III is very similar to the shape of the first mode 
in stage I, which is the first fundamental bending mode of the overhung rotor 
system. In addition, from figures 8 to 13 , it was found that none of the tor- 
sional vibratory frequencies and gear meshing frequencies are excited in the 
lateral rotor vibration of the system. 

Figure 14 depicts the gear mesh forces with effects of both imbalance and 
gearbox motion excitations. Note that the major component of these forces is 
the dc component at zero frequency which is due to the constant applied torque. 
The other sizable force component is at the tooth meshing frequency of 1800 Hz. 
The higher ac component which results in higher meshing frequency component 
between gears 1 and III, is because of the lower rotational stiffness in the 
stage III rotor. This can also be observed by the higher dc (zero frequency 
component) in stage III compared to stage II, as shown in figures 15 and 16. 

Only a very small force component can be observed at the gear mesh frequency 
(1800 Hz) and first mode frequencies of stage I (355 Hz), stage II (550 Hz), 
and stage III (280 Hz). Again figures 15 and 16 outline the modal rotational 
vibration characteristics of all rotor stages with both excitations. Note in 
figure 15 that other than the zero frequency component, the other sizable vibra- 
tion is at its own rotational modal frequency (i.e., 355 Hz at stage I, 550 Hz 
at stage II, and 280 Hz at stage III). Figure 16 depicts the second mode 
rotational vibration characteristics. Note that similar conclusions can be 
reached for the second modal frequencies of 1090, 1610, and 820 Hz, respec- 
tively, for stages I to III. 


SUMMARY OF RESULTS 

This paper presents a vibration analysis with the effects of gearbox 
motion and mass imbalance for multistage gear transmission. The analysis 
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combines gear mesh dynamics and structural modal analysis to study the trans- 
mission vibrations. The major content of this work can be summarized as 
fol lows : 

1. A comprehensive approach is developed to combine the nonlinear gear 
mesh dynamics with structural, lateral, and torsional vibration of the system 
to determine the global system response. 

2. The modal method transforms the equations of motion into modal coordi- 
nates to reduce the degrees of freedom of the system of equations. 

3. Gear force observations in both the time and frequency domains provide 
good insights into the source of dominating response forces. 

4. The magnitude of the gear mesh force is inversely proportional to the 
rotor stiffness of the geared system. 

5. The influence of the gearbox motion to system vibration is more pro- 
nounced in a stiffer rotor system. 

6. Gear tooth mesh frequency and torsional modal frequencies have substan 
tial effects on rotational but not on lateral vibrations of the system. 

7. Knowledge of modal amplifications under various excitations provide 
an understanding of the vibrational characteristics of the system. This knowl 
edge is crucial for designing transmissions with improved performance and 
durabi li ty. 
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TABLE I. - SYSTEM NATURAL FREQUENCIES 


Mode 

Stage I 

Stage II 

Stage III 

Torsional natural frequencies, Hz 

1 


MM 

MQM 

2 



HU 

Lateral natural frequencies, Hz 

1 

MM 

160 

110 

2 


189 

200 

2 

mm 

264 

160 
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Bearing 
supports - 7 ^ 
s' l\ 


Stage I 


Stage II 


kth stage 


Figure 1 —Coordinate system for gear mesh force. 


I . - | Stage III 

Figure 2 —Three-stage rotor-bearing gear system 


First lateral 

First torsional 


Second lateral 

Second torsional 



(a) Stage I. 



(a) Stage I. 



(b) Stage II. 



(b) Stage II. 



(c) Stage III. 

Figure 3. — First vibration mode of three- 
stage system. 


L I 

(c) Stage III. 

Figure 4. — Second vibration mode of three 
stage system. 
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(c) Case 2 in relation to time. 
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(e) Case 3 in relation to time. 
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(g) Case 4 in relation to time. (h) Case 4 in relation to frequency. 

Figure 8 — First lateral modal excitation on stage I. 
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Modal amplitude, in. Modal amplitude, in. Modal amplitude, in. Modal amplitude, in 
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(c) Case 2 in relation to time. 
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(d) Case 2 in relation to frequency. 
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(f) Case 3 in relation to frequency. 
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(g) Case 4 in relation to time. 


(h) Case 4 in relation to frequency. 


Figure 9. — First lateral modal excitation on stage II. 
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Time, ms 

(a) Case 1 in relation to time. 


(b) Case 1 in relation to frequency. 
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(c) Case 2 in relation to time. 
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(d) Case 2 in relation to frequency. 
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(e) Case 3 in relation to time. 
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(f) Case 3 in relation to frequency. 
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(g) Case 4 in relation to time. (h) Case 4 in relation to frequency. 

Figure 10 — First lateral modal excitation on stage III. 
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Modal amplitude, in. Modal amplitude, in. Modal amplitude, in. Modal amplitude, in. 
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(e) Case 3 in relation to time. 
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(g) Case 4 in relation to time. (h) Case 4 in relation to frequency. 

Figure 1 1 —Second lateral modal excitation on stage I. 
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Figure 12. — Second lateral modal excitation on stage II. 
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Figure 13. — Second lateral modal excitation on stage III. 
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Figure 15. — First torsional modal excitation on stage I. 
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Figure 16 — Second torsional modal excitation on stage I. 
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